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An unconventional superconducting phase is explored developing out of a non Fermi-liquid phase 
of the two-channel Anderson lattice model. It is characterized by a composite order parameter 
comprising of a local spin or orbital degree of freedom bound to triplet Cooper pairs with an 
isotropic and a nearest neighbour form factor. The superconducting transition temperature peaks 
in the moderate intermediate valence regime and vanishes at integral valence. The gap function is 
non analytic and odd in frequency, and a pseudo-gap develops in the conduction electron density of 
states which vanishes as \uj\ close to ui = 0. 
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Introduction. Heavy Fermion (HF) superconductivity 
has drawn much attention since the discovery of super- 
conductivity in CeCu2Si2 [|| which is likely characterized 
by an anisotropic order parameter with symmetry yet to 
be determined. It became apparent over the last decade 
that almost all HF materials are unstable with respect 
to magnetic or superconducting phase transitions, which 
either compete with each other or can even coexist as 
found in uranium based materials. Non-Fermi liquid be- 
haviour is often found in the vicinity of a so-called quan- 
tum critical point which is a point in parameter space 
where the ordering temperature is quenched to T = by 
a control parameter like pressure or doping ||. Doped 
CeCu6- x Ag x is one of the most prominent examples [Q. 

Two channel Anderson and Kondo lattice models, how- 
ever, also exhibit non-Fermi liquid behaviour in the para- 
magnetic phase driven by unquenched and fluctuating 
local degrees of freedom || ||. This phase is character- 
ized by a large residual resistivity and entropy, and ill 
defined electronic quasi-particles. Fermi liquid physics 
is restored by cooperative ordering or applied magnetic 
field or stress || Q . In particular, for a non Kramers dou- 
blet crystal field state in U 4+ or Pr 3+ ions, a two chan- 
nel Kondo effect is possible ||. The magnetic "spin" of 
the electrons serves as channel index in this case. UBC13 
and PrFePi2 are prominent candidates for a quadrupolar 
Kondo lattice description, as their enormous resistivity 
(> 100/if2cm) is removed only by phase transitions (su- 
perconductivity in UBei3, antiferroquadrupolar order in 
PrFePi2). Indeed, it has been shown that commensurate 
and incommensurate orbital ordering, as well as ferro- 
magnetism are possible in this model depending upon 
coupling strength and filling |Q . Evidence of a first order 
transition to an odd frequency pairing state in the Kondo 
limit (near integral valence) has been adduced, which is 
a singlet in both spin and channel indices, and no q value 
was preferred for the center of mass momentum (COMM) 
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In this paper, the possibility and nature of a super- 
conducting phase transition out of the NFL paramag- 



netic phase within dynamical mean field theory (DMFT) 
H [l(| will be explored. We show that in a intermediate 
valence (IV), two channel Anderson lattice a 2nd order 
transition to a triplet spin, triplet channel spin (StCt) 
order parameter with zero COMM will develop. The gap 
function is odd in frequency and singular, leading to a 
quasi-particle density of states that vanishes linearly with 
energy, in agreement with specific heat and spin lattice 
relaxation data for many HF materials. The order may 
arise form either two channel quadrupolar or magnetic 
ground states, and the transition temperature peaks in 
the moderate IV regime and vanishes at integral valence. 
The order parameter can be viewed as a composite of a 
bound spin and an even frequency triplet Cooper pair. 
Moreover, our analysis shows that isotropic and spatially 
extended states contribute to the superconducting phase, 
where the latter dominate in the IV regime. The binding 
of the fluctuating local degrees of freedom to the itiner- 
ant Cooper pair in the superconducting state is a natural 
consequence of the tendency to remove their residual en- 
tropy similar to a magnetic phase transition, a concept 
which was introduced first by Abrahams et. al. |H|| . 

Composite order parameter. The two-channel peri- 
odic Anderson model ^ 

A = E e * c L c w + E E * x i% + E 

+j2(-v 1+a v {cl^i a + h.c} (i) 
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describes the coupling of two degenerate conduction 
bands, labelled by a spin a and a channel index a = ±1 , 
to two localized doublets at each lattice site via hy- 
bridization matrix elements (— l) 1+a V. X are Hubbard 
operators, and e? the band dispersion. All energies 
will be given in units of A = V 2 nNf, the hybridiza- 
tion width. Since the model (|l]) conserves spin and or- 
bital quantum numbers, the symmetry of the conduc- 
tion electron pair operator is classified in four sectors, 
comprising of the product space of spin and channel sin- 
glets or triplets. Introducing the transposed bi-spinor 
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{k) — cg_|, Cj; ^J, the matrix elements 

of the tensor pair operator in the triplet /triplet (StCt) 
sector are given by 
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(2) 



where <x acts in the spin sector and f in the channel 
sector. The local nature of the DMFT permits only the 
Cooper pairs whose /c-dependent form factor S(k) trans- 
forms according to Ti , the trivial irreducible representa- 
tion. The pair expectation values (PA }, however, always 
must vanish due to Pauli's principle. As a consequence, 
the anomalous Green functions in this sector have to be 
odd in frequency. This peculiar kind of pairing was first 
suggested for He 3 Jl^] and later revitalized Q in the 
context of High-T c . Coleman et. al. Jl4| investigated the 
possibility of an odd-frequency pairing for a single chan- 
nel Kondo lattice model and found a staggered order pa- 
rameter with finite COMM Q. However, their mean-field 
solution breaks spin-rotational invariance. 

Abrahams et. al. noticed that the time derivative 
C 77 ' = (^:C 7 (t)c 7 ')| t _ is non-vanishing in the odd fre- 
quency superconducting phase and may possibly furnish 
an order parameter JO). The derivative ^:C 7 (r) is equiv- 
alent to the commutator [H, c 7 (r)] and consists of a prod- 
uct of two operators in the case of the two-channel pe- 
riodic Anderson model ([!]): a conduction electron pair 
operator and a local magnetic or channel spin operator. 
Thus, one arrives at a composite order parameter, which 
correlates local and itinerant degrees of freedom: forma- 
tion of band Cooper pairs is stimulated in the presence of 
a local pair resonating between definite spin and channel 
states with proper symmetry. Using equation of motion 
techniques on (|lj), we derive with the definitions 
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and - G c _ c _ , ,{iuj n ) is the Fourier transform of the 
anomalous Green function — (Tier (t)c r , ,) - 

\ ^ kaa ^ ' —ka'a'f 
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the exact relations for the symmetrized order parameter 
components for the StCt sector 



Oij = sign(A£) [( Si Pf) - (r 3 P^)] 



(5) 



= (\AE\Dij + sign(A£)T w ) /V 2 , (6) 

where i,j — x,y,z, i(j) indexes the spin(channel) 
space, P s,t = jj-J2k^(k)' l P T (k)icryiTyT'ip(—k) is the 
SsCt, and P*' 8 = ^Y,j:S{k)^ T (k)iT v i(T y a rf(-k) the 
StCs pair operator. Applying (0), the symmetrized pair 
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FIG. 1: (a) Pair susceptibility Xp C1 vs temperature in the 
triplet/triplet (StCt) sector for q = n (staggered) and q = 0. 
Parameters: AE — —3 and n c = 2. (b) Dependence of T c 
and characteristic temperature T* on n c , for AE — —A. 



expectation values Dij (Tij) are given by Dij(Tij) — 

YZa D aa>( T a£')[ i(T y <T i}crcr'[iTyT j } aa ,. Isotropic and 



modulated Cooper pairs with a form factor S(k) = 
contribute differently to the local composite order pa- 
rameter Oij in different regimes: in the IV regime, i. e. 
\AE = E a — E a \/A < I, modulated pairs dominate which 
is consistent with strong hybridization between /-shell 
electrons and conduction electrons, the isotropic pair 
contribution vanishing at the degeneracy point AE = 0. 
The composite order parameter Oij (|[) is invariant under 
exchange of the local doublets, and therefore supercon- 
ductivity can arise from either a magnetic or quadrupolar 
ground state. 

Pair Susceptibility. An extensive discussion of the 
physical nature of the possible order parameters in the 
StCt sectors is motivated by the observation, that the 
channel and spin exchanging local particle-particle (pp) 
/-Green function n^ c is purely odd in the incoming and 
outgoing frequencies for T — > 0, in addition to an l/\w\ 
singularity |L5| favouring an odd-frequency gap function. 
Eqn. (J^j connects the composite order parameter opera- 
tor to derivatives of the pair operators in frequency space 
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(iu) n )cryi(—iu] n ) (7 = kaa). Thus, 



the pair susceptibility contains the same information as 
the composite order parameter susceptibility. The irre- 
ducible pp interaction T u was obtained as a matrix in 
Matsubara frequency space 



^r" = [xZY 



[nf: 



Woe 



(7) 



by well established procedures within the DMFT j^, 
[j], [| [h]]. The local two particle band GF matrix 



Xi 



pp 



T,nX pp (q)/N is given by xf P c ( 







5n,mG c (iuj n )G c (— iuj n ), where G c denotes the local con- 
duction electron Green function, and Ilfc is the local 
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direct pp GF. Then, the pair susceptibility reads 



1 



1 - A, 



regular terms 



(8) 



where A m is the largest eigenvalue of the matrix M — 
\Jx pp { < l)jiT tt (fy\fx pp {. < i) m the odd frequency sector. 
c\m = En *w„Vx(9,0)(iw„)[|A)] n This DMFT-analog to 
the Eliashberg equation was first used by Jarrell and co- 
workers H . 

StCt Sector. The pair susceptibilities of all nine com- 
ponents of the tensor order parameter (||) in the StCt 
sector are equal, as expected on symmetry grounds. No 
superconducting instability was previously found for the 
two-channel Kondo-lattice model || since the model is 
restricted to integer rif valence. In the StCt however, 
T c is the largest for the IV regime, hence charge fluctua- 
tion driven rather the spin exchange induced. In Fig. [l]a, 
the pair susceptibility Xp is shown in the stable mo- 
ment regime for half-filling and for COMM \q\ = and 
\q\ = it. A sign change is found at T w 8 ■ 10~ 3 A which 
is an indicator for a second order phase transition in the 
triplet/triplet sector with an uniform order parameter. 

We investigated the pair susceptibility xf>* Ct m the IV, 
the Kondo and the stable moment regime for band-fillings 
between quarter and 3/4- filling and COMM of \q\ = 
and \q\ = ir. A superconducting transition was always 
obtained, but only for q — 0. T c is filling dependent for 
fixed coupling constant g = AffJ — A/|AI?| as depicted 
for the IV regime in Fig. [l]b. T c is also compared to the 
characteristic temperature of the lattice, T* as defined 
via a reduction of the effective moment /u^j(T) = T\{T) 
to 0.4 of its high temperature value. Hence, T* serves 
as a phenomenological Kondo lattice temperature (see 
also 0). The pair susceptibility in the singlet /singlet 
sector (SsCs) agrees with the one reported for the Kondo 
lattice model J6|. The instability, however, occurs always 
at lower temperatures compared to the one in the StCt 
sector. 

DMFT for the Superconducting Phase. The Nambu 
Green function is an 8 x 8 matrix in spin and orbital 
space. Assuming no directional coupling between spin 
and orbital degrees of freedom, the anomalous 4x4 self- 
energy matrix may be written as g(iuj ni k)a 2 T2 [n s d"ii c T] 
similar to He 3 p7[ , where n s and n c are constant unity 
vectors in spin and channel space , and the amplitude 
function g(iuj n ,k) has to be odd in frequency. This 
reduces the full problem to a standard 2x2 size for 
g(z, k) and the diagonal self-energy S. In order to de- 
rive DMFT equations with a purely local self-energy ma- 
trix S c (z), the anomalous self-energy is restricted to 
isotropic pairs, e. g. g(z,k) — g(z). G c (z) denotes the 
medium matrix in which the effective impurity is embed- 
ded. It is related to a generalized Anderson width ma- 
trix through A(z) = V 2 er 2 G c (z) <x 2 and has normal and 
anomalous components describing quasi-particle propa- 
gation and pair-creation and annihilation, respectively. 




FIG. 2: Spectral function of the quasi-particle (p e (to)) and of 
the anomalous Green function pf(lu), close to T c for AE = 
—2 and half filling in comparison with the spectral function 
slightly above T c (dotted line). The inset shows a blow up of 
the frequency interval \u>\ < 1CP 3 to demonstrate the p{ui) oc 
I a; I behaviour close to lu — 0. 



The band electron self-energy S c (z) = T 1 + GT 

is determined by the local T-matrix T. The diagonal 
elements of the T-matrix are given by the local quasi- 
particle scattering matrices Tn(z) = V 2 Gf(z), and ist 
anomalous contribution is calculated from 



T 12 (iu n ) = -V 4 



Uy(iuj n ,iuj m ]0)f(iuj m ) . (9) 



The negative sign takes into account that Ai2(z) = 
-V 2 f(z) where /» = [G c (z)] 12 and e(z) = [G c (z)} n . 
We obtained a solution for the superconducting phase 
close to T c , using the normal state self-energy T, e (z) — 
T e (z)/(1 + e(z)T e (z)), iterating only the anomalous T- 
matrix T s and medium /. The spectra of the quasi- 
particle and the anomalous Green functions are shown 
in Fig. ||. Since the the quasi-particle (qp) scattering 
rate T qp = ^smT,(—iS) tends to become strongly reduced 
close to lj = as expected, the full DMFT(NCA) turns 
out to be unstable in the SC phase which is related to 
the well known convergence problems in the Fermi liquid 
regime The resulting qp density of states p qp (oj) 

below T c is proportional to |a->| for small uj as depicted 
in the inset of Fig. |[ as can also be analytically shown 
using a temperature broadened Lorentzian as an ansatz 
for g(z) which reflects the 1/z singularity of H pp fill . 

Discussion and Conclusion. There is a controversial 
debate (see and reference therein) as to whether an 
odd- frequency solution with a COMM Q = is connected 
to a minimum or a maximum of the free energy for a 
second order phase transition and, hence, may be ther- 
modynamically unstable. Assuming a Fermi-liquid nor- 
mal phase and an infinitesimal odd-frequency anomalous 
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self-energy g(z) at T c , Heid showed that one indeed ob- 
tains an increase of the free energy in the superconduct- 
ing phase for Q = . This also holds for a non- Fermi 
liquid normal phase. 

Based on the self-consistent solution for the anoma- 
lous Green function, g(z) converges to a finite value such 
that always a gap is generated in the lattice Green func- 
tions. g(z) might be approximated in leading order by a 
Lorentzian with a finite weight A = y / L qp aT. The as- 
sumption of an infinitesimal value of g(z) for all frequen- 
cies is violated and hence Heid's theorem is not applicable 
in our case. Additionally, in the presence of an anomalous 
medium, the quasi-particle self-energy is also modified 
and, therefore, the feedback onto the effective site must 
be calculated self-consistently. We expect that the self- 
consistent solution for the one-particle Green function in 
the superconducting phase exhibits lim w _>o F qp — > 0. In 
this case, a non-analytic g{z) oc A/z, even with an small 
A, will be able to always produce a gap. Even though the 
DMFT(NCA) equations for the superconducting phase 
are numerically unstable, they indeed indicate a reduc- 
tion of the local free energy in the presence of a finite 
anomalous Green function and a tendency to reduce the 
scattering rate for the quasi-particles. Therefore, even 
though we cannot derive conclusively the type of the 
phase transition, we believe that theses arguments favour 
a second order transition in the StCt sector. This phase 
transition is associated with a minimum of the free en- 
ergy since an energy gap develops in the spectrum of the 
local band Green function indicating an energy gain by 
condensation. 
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FIG. 3: The phase diagram of the two-channel periodic An- 
derson model: the superconducting transition temperature 
T c , the antiferromagnetic transition temperature T m (AF), 
the spin-density wave transition temperature T m (SDW) are 
plotted versus the effective dimensionless coupling constant 
g = nA/\AE\. StCt superconductivity dominates the IV 
regime and the cross-over region to the Kondo regime. 

The phase diagram shown for n c — 2.2 in Fig. [| sum- 
marizes our investigation of superconductivity and mag- 
netism in the two-channel Anderson model. Supercon- 



ductivity dominates the IV regime and the corresponding 
order parameter has spin and channel spin triplet sym- 
metry. An instability in the SsCs reported in the two- 
channel Kondo lattice could be reproduced but occurs 
at much lower temperatures. Spin-density wave phase 
transitions take over in the Kondo-regime at g ss 1.3 cor- 
responding to AE = — 2.4A. The SDW wave-vector is 
continuously shifted towards nearest-neighbour antifer- 
romagnetism, which is suppressed for g — > 0. Since all 
calculations were performed in the paramagnetic phase 
of the model, the highest transition temperature defines 
the nature of the incipient order. One cannot rule out, 
however, that an SDW or orbital ordering phase is re- 
placed by, or even coexists with, a superconducting phase 
as found in some Uranium based HF compounds. We 
also would expect that the antiferroquadrupolar ordered 
PrFe4Pi2 might become superconducting when pressure 
is applied. 
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